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We study the influence of reheating on super-horizon density perturbations and gravitational 
waves. We correct wrong claimstl about the joining of perturbations at cosmological transitions 
and about the quantization of cosmological perturbations. 



1 Introduction 

The aim of these proceedings, based on Ref. 2, is to clear up the recent controversy on the 
relative contribution of density perturbations and gravitational waves to the CMBR today. 

The background model is taken to be a spatially flat FLRW model. We restrict our con- 
siderations to density perturbations and gravitational waves. The most general line element 
reads: 



ds^ = a^(r?){-(l + 20)dr?2 + 2B|idxM7? + [(1 - 2ij)jij + lE^nj + hj^dx'dx^} . (1) 
In the synchronous gauge, without loss of generality, the same line element can be written as: 

ds^ = a^(r/){-dr?2 + [(1 + hQ)6ij + ^Q,ij + /igwQij]dxMx^}. (2) 

The scalar function Q satisfies the Helmholtz equation and Qij is a symmetric, transverse and 
traceless spherical harmonic, k is the comoving wave number. In the following we use this 
form of the line element to make contact with previous works 0. It is convenient to define 
TC = a' /a, which is related with the Hubble constant H by the relation H = TC/a. A prime 
denotes a derivative with respect to conformal time. The quantity 7 is defined by the expression: 
7(7/) = 1- n'/n^. For the de Sitter space-time 7 vanishes. 

Let us now consider the equations of motion for the perturbed metric. For density pertur- 
bations (in the case of a vanishing anisotropic pressure) , all relevant quantities can be expressed 
in terms of the variable n defined by: ^ = a{h' -\- Tl^h)/{7{y/^), except in the de Sitter case, 
which must be treated separately. For gravitational waves, the relevant quantity is /Xgw = ahg^. 



Using the perturbed Einstein equations, one can show that both types of perturbations obey 
the same class of equation, i.e. the equation of a parametric oscillator: 

/i" + - = , (3) 

with Udp = (a-v/7)"/('^\/7) ^^'^ ^gw = a" /a. But there is a fundamental difference: the presence 
of the factor in the effective potential of density perturbations. 

Equation (^) is valid for any model. However, it is important to consider cases where exact 
analytical solutions can be found. This happens for power law inflation where the scale factor 
is given by a(r/) = lo\r]\^^^,P < —2. We will consider a model with three epochs in succession: 
inflation, radiation-dominated era, and matter-dominated era. For this simple model, the func- 
tion 7(r/) is a constant during each epoch and is given by: 7i = (2 + 71- = 2, 7m = 3/2. 
Therefore the constant factor ^/7 drops out of f/dp(^)) and the solutions for density perturbations 
and gravitational waves are given by the same Bessel functions: 

Kr]) = (A:r/)^/2[^i J/3+1/2 (A:??) + ^2 J-(;3+i/2) (^r?)]. (4) 

The aim is now to compute the amplitude of both types of perturbations during the matter- 
dominated era. In order to perform this calculation, two questions must be addressed: 

1) The initial conditions must be fixed. This amounts to choose the coefficients Ai and A2. 
This will be done with the help of quantum mechanical considerations. 

2) The way the solutions are matched between different epochs must be specified. This is no 
problem for gravitational waves since the effective potential is well-defined (although discontin- 
uous). This is more tricky for density perturbations. 7 is a Heaviside function. This means that 
the effective potential of the density perturbations is not defined in the sense of distributions 
at the different transitions. Recently, there was a controversy on this point §'i. One purpose 
of Ref. 2 and this paper is to clear up this question. This question arises because we consider 
simple models which allow analytical solutions. In reality, the transition is smooth and the 
effective potential is never ill-defined. 

2 Determination of the initial conditions 

Let us start with the first question. The normalization of the perturbed scalar field is fixed by 
the uncertainty principle of Quantum Mechanics. In the high frequency regime, this leads to: 

where the annihilation and creation operators satisfy the usual commutation relation. A hat 
indicates that we are now dealing with operators. The normalization of the scalar perturbations 
is fixed by the normalization of the perturbed scalar field since they are linked through Einstein's 
equations. In the high frequency limit, this link is expressed as: 

lim /i(?7, k) = — \/2Ka(^i(ry, k), (6) 

where k = SttG. For the model under considerations, it is easy to find the time dependence of 
the creation and annihilation operators. As a consequence all quantities are fixed uniquely, if 
at some initial time rjo the scalar field is placed in the vacuum state Ck(?y = %)|0) = 0. For 
gravitational waves the line of reasoning is similar. The resulfi is that the A's satisfy = 
I^J^^P^I = = l^^'^P)] = (27r/pi)/(|cos(7r/3)|\/^). We will show below that | A^'^^V^S'^^^ I = 1 

turns out to be crucial to predict the relative contribution of both types of perturbations. 



3 The junction conditions 



Let us now turn to the second question, i.e. the junction conditions. We analyze the approach 
taken by Grishchul0 and compare it with the claims of Deruelle and Mukhano^H. Suppose 
that the spatial transition hypersurface S is defined by its normal n^. To join two space-time 
manifolds along S without a surface layer the induced spatial metric hij = gij + n,nj and the 
extrinsic curvature Kij = (— l/2)£„/ijj should be continuous on S. In order to compute Ky the 
system of coordinates (i.e. the gauge) and the vector (i.e. the surface of transition) have 
to be specified. Different choices for lead to inequivalent junction conditions. If the surface 
of transition is defined by q{r],x^) = cte, it was shown by Deruelle and MukhanovH that the 
matching conditions read: 

[^ + n^h = o, [E]± = 0, (7) 



W + n^+{n' -n')-^]± = Q, [b-e' + -^]^ = q. (s) 
% % 

This result was generalized for spatially curved backgrounds and for the other types of perturba- 
tions in Ref. 2. If the surface of transition is a surface of constant time, the matching conditions 
in the synchronous gauge are: 

[h]± = [h'U = [hih = [h'i]± = 0. (9) 

For a sharp transition (i.e. if 7 is discontinuous) the two sets are not equivalent. Since it is 
believed that the transition occurs on a surface of constant energy which is not a surface of 
constant time, using the second set of joining conditions in this situation leads to an error. 
However, if the transition is smooth, i.e. \p]± = 0, then the two sets are equivalent, see Ref. 2. 

We use this result to study the regularisation scheme of Grishchulfl: Instead of matching 
inflation to the radiation epoch directly, Grishchuk introduced a smooth transition in between. 
Physically, this smooth transition represents the reheating of the Universe. It starts at 77 = 771 — e 
(end of inflation) and ends ai r] = rji + e (beginning of radiation). Here, e is small compared to 
r]i because we assume that reheating is a fast process. In the limit e goes to zero, we recover 
the sharp transition considered before and 7(7?) becomes a Heaviside function jumping from 
(2 -|- /?)/(! + P) to 2. For a smooth transition, without taking all the details of the reheating 
process into account, we do not know how the scale factor (and therefore 7) evolves between 
r]i — e and r/i -|- e. The idea of Ref. 1 was to assume that the function 7(7/) is given by: 

4 + 3/? , /3 ^ ^ ( r]-m \ i.r,. 

= WW) + WW) \-w) ' 

where s is a parameter controlling the sharpness of the transition. This equation holds for 
inflation and reheating, i.e. for rj between —00 and 771 + e. Formula ( |lO| ) leads to a reasonable 
equation of state p/ p = — 1 + 2/87(7/) and therefore gives a reasonable approximation to the real 
(exact) complicated function 7(77) even if details of the reheating process cannot be taken into 
account in such a simple approach. 

From the previous discussion, it is clear that 7(7/) is always continuous, even at r/ = 771 + e 
where the explicit joining was performed in Ref. 1 Ij. This means that \p\± vanishes. In Ref. 3, 
Deruelle and Mukhanov criticized the calculations done in Ref. 1 by means of the smooth 
transition described before, arguing that the junction conditions were not taken into account 
properly. We have shown that for continuous pressure the two sets of matching conditions are 

° In this paper rji was used instead of r;i + e to denote the end of reheating. It is very important to distinguish 
these two events. 



equivalent. Therefore the claim of Deruelle and Mukhanov is not appropriate. For a smooth 
transition, the matching conditions used by Grishchuk are perfectly justified since they coincide 
with the ones derived in Ref. 3. The argument of Deruelle and Mukhanov would be relevant if 
the transition was sharp and 7(7/) discontinuous at rj = rji + e. 

Even the simple form (|l^) is too complicated to allow a direct integration of the equation 
of motion for //. Nevertheless, we can follow the evolution of // through inflation and reheating. 
For rj < rji — e, 7(r/) is a constant and the solution for fj, is given by Eq. (Q). The value of fi{r]) 
just before reheating can be easily computed: 

4 (dp) 4 (dp) 



2^+2r(/3 + |) 2^+2r(/? + |) 

because krji <C 1 and e <^ r]i. Between rji — e and r/i + e the function 7(7/) is no longer a constant 
and the solution (^) can no longer be used. In order to evolve /i through the reheating transition 
we use the superhorizon solution of Eq. @, /U ^ 0-^7, to obtain 




^(77i + e)~ — . , . a(?7i +6)^7(7?! + g) ~^(7/i-e)J- , (12) 

- e)v 7(^/1 - e) V 7i 

because a(ryi + e) a(r/i — e). This relation should be compared to Eq. (81) and to the 
relation /x|^^_o = /"|r?i+o below Eq. (48) of Ref. 1. From Eq. (|T^), it is clear that the ratio 
/j,{rji + e)//i(??i — e) is not 1 but proportional to 1/^/7^. This factor is huge when 7i is close to 
(de Sitter). Therefore the mistake in Ref. 1 was not the use of wrong junction conditions but the 
fact that ^(77) was not evolved correctly through the reheating transition: 7(771 — e) 7^ 7(^1 + c) 
implies iJ.{rii - e) / /i(r/i + e). 

Let us turn to the radiation-matter transition taking place at equality. In Ref. 1 the equality 
transition was treated as a sharp transition. The joining conditions (^) are not correct in 
general (i.e., for any residual gauge fixing). However, if one specifies the synchronous gauge in 
the matter dominated epoch to be the comoving one, then the joining conditions are fine at 
the equality transitional. That is because the density contrast vanishes at the leading order [i.e., 
it is proportional to (A;r/)^]. 

Using the previous results, it is interesting to calculate the ratio of the Bardeen potential 
<!> = (/>+ [a{B — E')\' /a and /igw at superhorizon scales today. We obtain (see Ref. 2): 

10V2A(^")M^i-6) _ 10 ^^^^ 



today 3 AfP)/^(m+e) 3 



This equation illustrates the importance of the initial conditions and the importance of the 
evolution of density perturbations during reheating. For small values of 7i the amplitude of 
scalar metric perturbations is larger than the amplitude of gravitational waves, which implies 
that the main contribution to the large scale fiuctuations in the cosmic microwave background 
is due to density fluctuations. 
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Titre en frangais: Amplitude des perturbations cosmologiques plus grandes que I'horizon. 



Resume en frangais: 

Nous etudions I'influence du 'reheating' sur les perturbations de densite et les ondes gravi- 
tationnelles plus grandes que I'horizon. Nous corrigeons de fausses affirmations faites a propos 
du raccordement des perturbations lors des transitions cosmologiques et a propos de leur quan- 
tification. 



